Nuclear target effect on dark matter detection rate 
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Expected event rates for a number of dark matter nuclear targets were calculated 
in the effective low-energy minimal super symmetric standard model, provided the 
lightest neutralino is the dark matter Weakly Interacting Massive Particle (WIMP). 
These calculations allow direct comparison of sensitivities of different dark matter 
detectors to intermediate mass WIMPs expected from the measurements of the DArk 
MAtter (DAMA) experiment. 

I. INTRODUCTION 

Weakly Interacting Massive Particles (WIMPs) are among the most popular candidates for 
the relic cold dark matter (DM). There is some revival of interest to the WIMP-nucleus spin- 
dependent interaction from both theoretical (see e.g. and experimental 



sec c. 



18L ll9L l20j ) points of view. There are some proposals 



aimed at direct DM detection with relatively low-mass isotope targets 



as well as some first attem pts to desi gn a nd construct a DM detector which is sensitive to 



the nuclear recoil direction 
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26, 



271 ] . Kinematically, low- mass targets make 
preference for the low-mass WIMPs (due to Mrarget ~ Mwtmp) and are more sensitive to 
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13- On 



the other 



the spin-dependent WIMP-nucleus interaction as well jl, 
hand, WIMPs with masses about 100 GeV/c 2 follow from the results of the DArk MAtter 
(DAMA) experiment. This collaboration claimed observation of the first evidence for the dark 
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Aimed for 



matter signal due to registration of the annual modulation effect |3lL 
more than one decade at the DM particle direct detection, the DAMA experiment with 100- 
kg highly radio-pure Nal(Tl) scintillator detectors successfully operated till July 2002 at the 
Gran Sasso National Laboratory of the I. N.F.N. On the basis of the results obtained for over 7 
annual cycles (107731 kg-day total exposure) the effectiveness of the WIMP model-independent 
annual modulation signature was demonstrated and the WIMP presence in the galactic halo is 
strongly supported at 6.3 a C.L. j^. The main result of the DAMA observation of the annual 
modulation signature is the low-mass region of the WIMPs (40 < M WIMP < 150 GeV/c 2 ), 
provided these WIMPs are cold dark matter particles. It is obvious that such a serious claim 
should be verified by other independent measurements. To reliably confirm or reject the DAMA 



2 



result an experiment should have the same or better sensitivity to the annual modulation signal. 
At the same time one should know an expected rate in his own detector, provided the DAMA 
result is correct. Motivated by the DAMA evidence predictions for the direct DM detection 
rate in a Ge-73 detector within the framework of the so-called effective low-energy minimal 
super symmetric standard model (effMSSM) are given in Here the analysis of is 

extended and the expected WIMP detection rates are recalculated in the effMSSM for other 
DM-interesting isotope targets on the basis of the data base from j^J . Similar expectations are 
obtained as well by means of direct recalculations of the DAMA constraints j^] into detection 
rates for other targets. This allows one to estimate the prospects to confirm or to reject the 
DAMA result by other DM search experiments. 



II. EVENT RATES AND EFFECTIVE LOW-ENERGY MSSM 



A dark matter event is an elastic scattering of a relic WIMP \ from a target nucleus A 
producing a nuclear recoil E R which can be detected by a suitable detector. The differential 
event rate in respect to the recoil energy is the subject of experimental measurements. The 
rate depends on the distribution of the WIMPs in the solar vicinity f (y) and the cross section 
of WIMP-nucleus elastic scattering l4lH . The differential event 

rate per unit mass of the target material can be given as a sum of spin- dependent (SD) and spin- 
independent (SI) contributions, parameterized via spin-independent (cr SI ) and spin-dependent 



<7 SD ) WIMP-nucleon interaction cross sections 



32, 



3J: 



dR(E R ) A7 p x /""max do 2 , . , . 

- N T — / dvf{v)v—{v,q ) = K sl (E R ,m x ) a S i + K SD (E R ,m x ) a SD . 



dE R m x J Vmin dq 1 

The nuclear recoil energy E R = g 2 /(2Myi) is typically about 10 _6 m x . The number of nuclei 
per unit of target mass is Nt and Ma is the target nucleus mass. The effective spin WIMP- 
nucleon cross section cr SD and the coupling mixing angle 9 were introduced 
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421 in such a 



way that SD WIMP-proton and SD WIMP-neutron interaction cross sections have the form 



a, 



SD 



cr SD • cos 2 9, and aUr, = cr QT ^ ■ sin 2 9. Further notations are: 



'SD 



'SD 



K sl (E Rl m x ) = N T ^^B sl (E R ) 
zm x [i p 

Nt ^ — ^ B SD (E R ) 



Ml 



K SD {E R ,m x ) 
Bsi,sd(E r ) 



2m x nl 



4 J + l 
3 J 



((Sp) cos 9 + (S„) sin^y 



(2) 



(v 2 ) 



Here (S p ( n )) is the spin of the proton (neutron) averaged over all nucleons in the nucleus A. 
The dimensionless integral I(E R ) is dark-matter-particle velocity distribution correction, which 
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reduces the rate for large enough momentum transfer: 

r,^ (v 2 ) f f( x ), V* 3 + 2ry 2 . . 

I{Er) = tt / dx = — — — a erf(g miw + 77 - erf x min - 77 , 

W Jimi, v 2 ^(1 + 2r/ 2 )erf(?7) + 2rje *> 

where we assume that in our Galaxy rest frame WIMPs have the Maxwell-Boltzmann velocity 

distribution and use the dimensionless Earth speed with respect to the halo 77 = 1, x 2 mia = 

— — f„ R , ua = — 7~r- The error function is erf fx) = —= [ dte~ t2 . The velocity variable 

4 n A v 1 m x + M A y/n Jo 

is the dispersion v ~ 270km/c. The mean WIMP velocity (v) = J\v. For the WIMP mass 
density in our Galaxy the value p x = 0.3 GeV/cm 3 is used. We also assume both form-factors 
Fgj sj)(Er) in the simplest Gaussian form following 0, 44]. In particular, this allows rather 
simple formulas (0) to be used, which are suitable for our comparative consideration. The total 
direct detection rate R(e, e) can be obtained by integrating differential rate over the recoil 
energy interval from the threshold energy e to the maximal energy e. To accurately estimate 
the event rate R(e, e) one needs to know a number of quite uncertain astrophysical and nuclear 



structure parameters as we 
for example, discussions in 



1 as the very specific characteristics of an experimental setup (see, 
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45(). In this paper it is enough to assume these uncertainties to 



be almost the same for all target materials considered. Furthermore, one should calculate cross 
sections cr gD and er gi within the framework of, for example, some SUSY-based theory or take 
them from experimental data. In both cr SD and <7 gi have already been calculated within 
a phenomenological SUSY model whose parameters are defined directly at the electroweak 
scale (see e.g. This effective scheme of 

the MSSM is called the effMSSM in 35], and later the low-energy effective supersymmetric 
theory (LEEST) in (5^, 0|. The effMSSM parameter space is determined by entries of the 
mass matrices of neutralinos, charginos, Higgs bosons, sleptons and squarks. The list of free 
parameters includes tan (3, the ratio of neutral Higgs boson vacuum expectation values; /i, the 
bilinear Higgs parameter of the superpotential; M lj2 , soft gaugino masses; M A , the CP-odd 
Higgs mass; m~, m~, m~ (m~, m~), squark (slepton) mass parameters squared for the 1st and 
2nd generation; m~ , m~, m~ (m~ , m~), squark (slepton) mass parameters squared for the 
3rd generation; A t , A&, A T , soft trilinear couplings for the 3rd generation. The third gaugino 
mass parameter M 3 defines the mass of the gluino in the model and is determined by means of 
the GUT assumption M2 = 0.3 M3. The intervals of the randomly scanned MSSM parameter 
space in j3| were narrowed to fit the DAMA-inspired domain of the lower masses of the LSP 
(m x < 200 GeV). The current experimental upper limits on sparticle and Higgs masses from 
the Particle Data Group as well as the limits on the rare b — >• 37 decay have been imposed. 
For each point in the MSSM parameter space (MSSM model) the relic density of the light 



neutralinos ^l x h^ was evaluated with the code 



Q,Q,Q 



based on DarkSUSY 



56 with the 



allowance for all coannihilation channels with two-body final states that can occur between 
neutralinos, charginos, sleptons, stops and sbottoms as l ong as their masses are rrii < 2m x . 
Two cosmologically interesting regions are considered in |34j ]. One is 0.1 < Q x h 2 < 0.3 and 



the other is the WMAP-inspired region 0.094 < Q x h 2 < 0.129 57, 58]. Further details can be 



found in 



[34]. 



0Q. 



III. RESULTS AND DISCUSSION 
A. Calculations in effMSSM 



Integrating the differential rate from the recoil energy threshold e to some maximal 
energy e one obtains the total detection rate R(e, e) as a sum of the SD and SI terms: 

R(e, e) = R S i(e, e) + R SD (e, e) = J dE R K sl (E R , m x ) a S i + J dE R K SD (E R , m x ) a sr> . (3) 



In |3J estimations of the ideal total expected rate i?(0,inf) for WIMPs with M WIMP < 
200 GeV/c 2 in a 73 Ge detector are obtained within the effMSSM. Here, with cx SD and o~ SI 
already calculated in j^j], new estimates of the integrated event rate R(e,e) for WIMP masses 
smaller than 200 GeV/c 2 are obtained for a number of DM targets. For definiteness, the recoil 
energy threshold e = 5 KeV (and sometimes e = 10 KeV) with the maximal energy e = 50 keV 
are used. 

The calculated event rates _R(5, 50) and the rate ratios for different targets are depicted as 
scatter plots in Figs. HHS1 For example, in Fig. [T]one can see total, i?(5,50), spin- independent, 
R(5, 50)si and spin-dependent _R(5, 50)sd event rates expected in 73 Ge-target together with 
their ratio -R(5, 50)sd/-R(5, 50)si as functions of the WIMP mass. It is always interesting to 
trace some interplay between SD and SI contributions to the total event rate. To this end 
correlations between the total rate -R(5, 50) and its SD fraction as well as correlations between 
the total rate i?(5, 50) and the SD-to-SI ratio _R(5, 50)sd/-R(5, 50)si are also given in Fig. [T] All 
open symbols in the figures correspond to the case when the WMAP constraint on the relic 
neutralino density 0.094 < Q x h 2 < 0.129 is taken into account. The relevant filled symbols show 
the rates which one would expect if an extra DAMA constraint, 1 ■ 10~ 7 pb < Ogj < 3 ■ 10 -5 pb, 
is imposed on SI cross sections. It is seen that in the last case (filled symbols) the SI rates 
are at least two orders of magnitude larger than the SD one and the large total rate values 
(R > 0.01 events/day/kg) are saturated only by the SI interactions. If one ignores these filled 
symbols (i.e. the DAMA-inspired extra constraint 1 • 10 -7 pb < Ogj < 3 • 10~ 5 pb), then the 
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FIG. 1: Expected in 73 Ge total i?(5, 50), SI and SD event rates i?(5, 50)si and R(5, 50)sd (upper panel) 
as well as the ratio R(5, 50)sd/-R(5, 50)si (lower left panel) as functions of WIMP mass. Correlations 
between the total rate R(5, 50) and the SD fraction in R(5, 50) is given in the lower middle panel. 
Correlations between the total rate R(5,50) and the ratio R(5, 50)sd/-R(5, 50)si is given in the lower 
right panel. Open symbols correspond to the WMAP constraint 0.094 < £l x h? < 0.129. Closed 
symbols give rates with an extra DAMA constraint 1 • 10~ 7 pb < erigj < 3 • 10~ 5 pb. 
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SD contribution does not look very suppressed and the SD contribution alone can saturate the 
total event rate, but only when the rate itself is rather small (R ~ 0.001 events/day/kg). These 
features take place for all heavy enough targets, therefore the corresponding figures for Nal, 
Csl, and Xe target are not given. 

It is well known that a fluorine-containing target is the best one for detection and measure- 
ment of the spin-dependent WIMP-nucleus interaction (see e.g. Figure El shows that 
the SD rate in CF 4 is indeed the biggest one and for a large number of points (MSSM models) 
the SD contribution dominates. Nevertheless it is also seen that it is not correct to completely 
ignore the SI contribution to the total expected rate in the fluorine target U [y| because the 
SI rate is almost the same as the SD one. Furthermore, at a current level of the DM detector 
sensitivity, when the DAMA-inspired large SI contributions are not yet completely excluded 
(filled symbols in Fig. EJ), the SD contribution in CF 4 , i?(5,50)sD, is smaller than the SI one, 
i?(5, 50)si. The ratios of the total, SI and SD rates in the CF 4 and 73 Ge targets are presented 
in Fig. El function of WIMP mass. 




WIMP mass, GeV 



50 200 
WIMP mass, GeV 



FIG. 3: Ratios R(5, 50)cf 4 /-R(5, 50)73Q e of the total (top panel), SI (middle) and SD (bottom) rate 
as functions of WIMP mass. No relic density constraint is imposed in the left column. Open symbols 
in the right column correspond to the WMAP constraint 0.094 < Q x h? < 0.129. Filled symbols give 
these ratios for the rates obtained with an extra DAMA SI constraint, 1 • 10 -7 pb < crgj < 3 • 10 -5 pb. 



The ratios in the left panel of Fig.Elcorrespond to the rates calculated without any constraint 
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FIG. 4: Expected in the CH4 target the total R(5, 50), SI and SD event rates i?(5, 50)si and R(5, 50)sd 
(upper panel) as well as the ratio R(5, 50)sd/-R(5, 50)si (lower left panel) as functions of WIMP mass. 
Correlations between the total rate i?(5, 50) and the SD fraction in i?(5, 50) is given in the lower 
middle panel. Correlations between the total rate R(5, 50) and the ratio i?(5, 50)sd/-R(5, 50)si is given 
in the lower right panel. Open symbols correspond to the WMAP constraint 0.094 < Q x h 2 < 0.129. 
Closed symbols give rates with an extra DAMA SI constraint, 1 • 10~ 7 pb < cigj < 3 • 10 -5 pb. 
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on the relic density of neutralinos in the effMSSM. The increase of these ratios at very low WIMP 
masses reflects better sensitivity of fluorine to smaller WIMP masses than that of germanium. 
Open symbols in the right panel of Fig. El depict these ratios when the WMAP constraint 
0.094 < Q x h 2 < 0.129 is imposed on the calculated neutralino relic density. Filled symbols 
gives these ratios for the rates obtained with an extra DAMA SI constraint, 1 ■ 10~ 7 pb < Ogj < 
3 • 10~ 5 pb. The sensitivity of CF 4 to the SD WIMP interaction is about ten times as large as 
the SD sensitivity of 73 Ge. At the same time, the CF 4 sensitivity to the SI WIMP interaction 
is less than 0.1-0.05 of the SI sensitivity of 73 Ge. As a result, the total expected rate in 73 Ge is 
a bit larger than in a very spin-sensitive CF4 target. The right panel in Fig. El shows that the 
relic density WMAP and extra DAMA constraints make this conclusion stricter. The expected 
total rate in a heavy enough 73 Ge target is about ten times as large as the total expected rate 
in the CF 4 target. 




50 100 150 200 50 100 150 200 

WIMP mass, GeV WIMP mass, GeV 



FIG. 6: Ratios R(5, 50)3 He /-R(5, 50)73 Ge of the total (top panel), SI (middle) and SD (bottom) rate as 
functions of WIMP mass. No relic density constraint is imposed in the left column. Open symbols 
in the right column correspond to the WMAP constraint 0.094 < Q x h? < 0.129. Filled symbols give 
these ratios for rates obtained with an extra DAMA SI constraint, 1 • 10~ 7 pb < <7gj < 3 • 10~ 5 pb. 

Figures El and El show the total, SD and SI expected event rates in the CH4 and 3 He targets 
which contain the lightest nonzero-spin nuclei, interesting for direct DM search. Absolute 
values of the total, SD, and SI rates in these materials are visibly smaller than in the fluorine- 
containing CF4 and germanium targets (Figs. EHZj); especially when the extra WMAP (left 
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50 100 150 200 50 100 150 200 



WIMP mass, GeV WIMP mass, GeV 

FIG. 7: Ratios R(5, 50)ch 4 /-R(5, 50)73g c of the total (top panel), SI (middle) and SD (bottom) rate 
as functions of WIMP mass. No relic density constraint is imposed in the left column. Open symbols 
in the right column correspond to the WMAP constraint 0.094 < Q x h 2 < 0.129. Filled symbols give 
these ratios for rates obtained with an extra DAMA SI constraint, 1 • 10 -7 pb < Ogj < 3 • 10~ 5 pb. 
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FIG. 8: The same as in Fig. EJ but for Nal and 73 Ge. 



10 



panel) and DAMA (filled symbols) constraints are imposed. 

For all detectors with heavy enough target mass (Nal, Csl, Xe, etc) the absolute values of 
the total, SI and SD rates look very similar with the only possible exception in the domain 
of very low-mass WIMPs (when a target contains some light isotope like, for example, Na in 
the Nal target). In Figure BDa set of NaI-to- 73 Ge rate ratios, R(5, 50)^/ R(5, 50)73 Ge , is given 
for illustration of the behavior. All depicted ratios are of the order of unity. Only for very 
low WIMP masses (less than 10 GeV/c 2 ) and with the WMAP constraint neglected the rates 
in Nal start to clearly dominate over the rates in 73 Ge due to kinematically preferable WIMP 
interaction with Na. This low-WIMP-mass growth of the ratios is absent in the Csl and Xe 
targets. Furthermore, for these materials the rates in 73 Ge dominate in the low-mass WIMP 
region. 



B. Calculations with DAMA constraints 



In the previous section the WIMP-nucleon a SD and a si cross sections which enter event rate 
(JHJ) were taken from theoretical calculations in the effMSSM [34]. Another source of these cross 
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FIG. 9: Variations of expected event rates, i?(5, 50), for a number of targets followed from the DAMA- 
allowed cross sections c gD and <t si . Targets with nonzero-spin nuclei from the odd- neutron (odd- 
proton) group model are given in the left (right) part of the figure. 



sections is, for example, the DAMA experiment 32]. With formulas (£Q)-(HD) the DAMA(Nal) 
constraints on o" SD and <r SI from j^] can be phenomenologically transformed into the allowed 
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FIG. 10: Variations of expected spin-independent contributions to the event rate, R(5, 50)si, in a 
number of targets followed from the DAMA- allowed cross sections crg D and <7 SI . 



T3 

& 10 ? 



B 1 
> 



o 



a 10 

a 

10 



10 



1 



10 



'He 



Ge-nat 



Ge 



"Xe 



"Xe 



CF, 



GaF, 



ALO, 
2 Nal. 



FIG. 11: The same as in Fig. 1101 but for the spin-dependent contributions -R(5,50)sd- 

regions for the detection rates with other targets. The results of these recalculations are given 
in Figs. IUHTT1 for the total, SI and SD expected rates R(10, 50) in a number of representative 
materials for a DM detector. Here the threshold of 10 keV is used. The values of the expected 



12 



rate can vary within the columns without any conflict with the DAMA-allowed <r SD and cr gI 
regions. The left (right) parts of these figures contain rate restrictions for the odd-neutron 
(odd-proton) group model nonzero-spin nuclei (see e.g. 0). Some of the highest rate values, 
for example in 73 Ge, are already excluded by measurements In particular, from Figs. IUHTT1 
one can see that all fluorine-containing targets (LiF, CF4, C2F6 and CaF2, etc) have almost 
the same sensitivity to the both SD and SI WIMP-nucleus interactions. Among all materials 
considered a detector with a 73 Ge, 129 Xe, or Nal target has better prospects to confirm or 
to reject the DAMA result due to the largest values of the lower bounds for the total rate 
(i?(10, 50) > 0.06 — 0.08 events/day/kg). If, for example, one ignores the SI WIMP interaction 
(Fig. HOj) . then all materials have almost the same prospects to detect DM particles with the 
only exception of CH 4 . 



IV. CONCLUSION 



Expected event rates for a number of dark matter target materials are calculated in the 
effective low-energy minimal supersymmetric standard model (effMSSM), provided the light- 
est neutralino is the dark matter Weakly Interacting Massive Particle (WIMP). The results 
obtained are based on previous evaluations of the neutralino-proton (neutron) spin and scalar 
cross sections for the neutralino masses m x < 200 GeV/c 2 [34]. The performed calculations 
allow direct comparison of sensitivities of different dark matter setups to the WIMPs expected 
from the measurements of the DAMA experiment. In particular, it is shown that detectors 
with a 73 Ge, 129 Xe, and Nal target have better prospects to confirm or to reject the DAMA 
result. 

It is worth noting, that to get very accurate predictions for the event rate one has to take 
into account a number of quite uncertain astrophysical and nuclear parameters and specific 
features of a real setup. We considered only a simple spherically symmetric isothermal WIMP 
velocity distribution 0, |^ and do not go into detail of any possible and in principle im por tant 
uncertainties (and/or modulation effects) of the Galactic halo WIMP distribution [21I l6l[ lo^l 



63, 



ft 
|67 



65, 66]. For simplicity we use the Gaussian scalar and spin nuclear form-factors from 



44l 1671 ] . We believe it is relevant for our comparative study because the very influence of the 



factors is suppressed in the rate ratios. 

The authors thank the RFBR (02-02-04009) and the JINR-Slovak Republic program for 
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